function y = initial_function1(x)
y = sin(x) + sin(10*x);
    